Generation of J- wave coupling in the two-dimensional Hubbard 
model from functional renormalization 
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Within the two-dimensional repulsive f- /'-Hubbard model, an attractive coupling in the li-wave pairing chan- 
nel is induced by antiferromagnetic fluctuations. We investigate this coupling using functional renormalization 
group equations. The momentum dependent rf-wave coupling can be bosonized by the use of scale dependent 
field transformations. We propose an effective coarse grained model for the Hubbard model which is based on 
the exchange of antiferromagnetic and rf-wave collective bosons. 
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I. INTRODUCTION 

The Hubbard model iH Q, Htl for strongly correlated elec- 
trons has been used for a wide variety of phenomena, ranging 
from high temperature superconductivity |4] to the metal in- 
sulator transition or antifeiTomagnetism. Solving this model 
is a major theoretical challenge. This is due to the com- 
plexity of the effective electron interaction, typically char- 
acterized by the competition of different channels, like the 
antiferromagnetic or (/-wave Cooper-pair exchange channels 
lH II ill 113, [H [12 [IS . So far, the development of 
d-wa\e superconductivity is still a controversial issue. Al- 
though several many-body techniques predict the emergence 
of a d-wave instability as leading instability in certain param- 
eter ranges, see e. g. lUt], numerical studies have difficulties 
to detect superconductivity ||l|ll|ll[ll[Il[Ill'l9',20], as a 
consequence of both finite size and temperature limitations. 

The c/-wave-pairing interactions are absent in the micro- 
scopic Hubbard model. As a plausible mechanism lead- 
ing to d-wave pairing the exchange of effective antiferro- 
magnetic bosonic degrees of freedom has been proposed 
iH 0, 1 H [H, [H. This idea was also suggested by Hlfll 
on the basis of a phenomenological spin-spin susceptibility. 
The generation of the interaction in the cZ-wave channel has 
been investigated extensively by the functional renormaliza- 
tion groiip study of a momentum dep endent four electron vst- 
tex iB E III IM IM 123 IM IM IM El. our inves- 
tigation employs non-perturbative flow equations, based on 
exact renormalization group equations for the average action 
or flowing action |32D, 113 3fl . It is complementary to studies 
of the electron vertex by focussing explicitely on the role of 
the antiferromagnetic bosons, using the technique of partial 
bosonization during the flow |[34l l35n . We investigate the two- 
dimensional Hubbard model with nearest and next to nearest 
neighbor hopping. 

The aim of this note is the derivation of an effective coarse 
grained model, valid at intermediate length scales fc"' larger 
than the lattice distance a, but smaller than the scale of macro- 
scopic physics. Our effective model will be based on a de- 
scription of the electron interactions by the exchange of col- 
lective bosons. Near half-filling a model with exchange of an- 
tiferromagnetic boson has been used ll36[|37ll for a quantitative 



computation of effective antiferromagnetic order below the ef- 
fective critical temperature Tc- It is advocated that in two di- 
mensions the size of the ordered domains increases exponen- 
tially as the temperature decreases towards zero. For T < Tc 
this size exceeds the size of a typical experimental probe. For 
all practical purposes the physics is then the same as for an or- 
dered system. This includes the existence of long range flucta- 
tions associated to Goldstone bosons. For Tc < T < Tpc short 
range antifeiTomagnetic order is found, with a typical domain 
size smaller than the experimental probe. Here Tpc denotes the 
pseudo-critical temperature below which short range antifer- 
romagnetic order sets in. In this paper we extend this approach 
by including the exchange of collective bosons consisting of 
Cooper-pairs in the d-wa\e channel. They are assumed to be a 
crucial ingredient for the understanding of the Hubbard model 
away from half filling. Their condensate would lead to super- 
conductivity. 

We have investigated earlier the phase transition to super- 
conductivity in an effective li-wave exchange model fsS]. The 
pha se transition is found to be of the Kosterlitz-Thouless type 
03911 . characterized by (modified) essential scaling above Tc, a 
jump of the superfluid density at Tc, and a gapless excitation 
with temperature dependent anomalous dimension below Tc- 
The present note constitutes a step for establishing such an ef- 
fective model as a coarse grained version of the microscopic 
Hubbard model. This mapping is not yet complete, since we 
concentrate here only on the generation of the effective elec- 
tron interaction in the li-wave channel and its bosonization, 
while a more detailed investigation of the propagators and in- 
teractions of the d-wa\e bosons is postponed to a subsequent 
publication. 

The partially bosonized model opens the door for a straight- 
forward computation of the flow below the pseudo-critical 
temperature Tpc, where local order sets in and higher order 
boson interactions (e.g. corresponding to eight-fermion inter- 
actions) play an important role. The low temperature region is 
notoriously difficult for a purely fermionic description since 
the four fermion interaction diverges as the pseudo-critical 
temperature Tpc is approached from above. Also the impor- 
tant effect of higher fermionic operators can be taken into ac- 
count most conveniently in a partially bosonized formulation 
ll36[|37ll . In the partially bosonized description the divergence 
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of the four fermion interaction is due to a boson "mass term" 
or "gap" changing from positive to negative values during the 
flow. A negative mass term indicates local order, since at a 
given coarse graining scale k the eff'ective potential has a min- 
imum for a nonzero value of the boson field. If this order per- 
sists for k reaching a macroscopic scale, the model exhibits 
effectively spontaneous symmetry breaking, associated in our 
model to antiferromagnetism or li-wave superconductivity. 

On the other hand, the partially bosonized formulation in- 
troduces some bias into the renormalization group procedure 
by focussing on one or several boson exchange channels. This 
can be circumvented in the purely fermionic renormalization 
group flow. However, this bias can be systematically reduced 
by the introduction of further bosonic fields, which yield addi- 
tional resolution for the description of the interaction between 
the electrons. The present approach constitutes a first step 
into this direction and is thus complementary to references 

il 111 111 IM 111 111, 111 111 111 [B III . 

II. EFFECTIVE AVERAGE ACTION 

We use a functional integral framework and investigate the 
effective average action SIH], which includes all quan- 
tum and thermal fluctuations with momenta > k^. This is 
realized by introducing an infrared cutoff which suppresses 
the fluctuations with < k^. This cutoff will be removed at 
the end for A: — > 0. For ^ = one recovers the usual ef- 
fective action - the generating functional for the one particle 
irreducible Green functions. For A: — > oo, or A; — > A, with A 
a suitable ultraviolet scale, the fluctuation effects are negligi- 
ble and Fa becomes the microscopic or classical action S a- 
Therefore the scale dependent average action F^ interpolates 
between the microscopic action 5 a and the full quantum ef- 
fective action F, 

limF^:^5A, limFt^F. (1) 

it— >A k—>0 

The A;-dependent flow of the average action obeys an exact 
renormalization group equation ll32ll . 

dkTkM = ^STr{[rf M + RkV'dM , (2) 

where the "supertrace" STr runs over field type, momentum 
and internal indices, and has an additional minus sign for 
fermionic entries. The functional differential equation ^ in- 

(2) 

volves the full inverse propagator F^ [x] (second functional 
derivative of F^), regulated by 7?^. Approximations to the so- 
lution of Eq. (|2]l proceed by a truncation on the rhs with a 
suitable ansatz for the form of F^. 

On the microscopic level, the Hubbard model is defined as 
a purely fermionic model for electrons on a lattice. How- 
ever, it is well known that also bosonic degrees of freedom, 
e.g. Cooper-pairs, play an important role on larger length 
scales. Hence, the question arises how the optimal descrip- 
tion of the relevant degrees of freedom on all scales can be 
achieved within a functional renormalization group study. On 
the one hand, we want to accommodate the purely fermionic 



model at high momentum scales. On the other hand, we want 
to study phase transitions and critical behavior for the macro- 
scopic physics. In the vicinity of the critical temperature the 
long range fluctuations are dominated by bosonic compos- 
ite operators. In order to exploit dkectly the importance of 
the collective bosonic degrees of freedom, we employ partial 
bosonization as motivated by a Hubbard-Stratonovich trans- 
formation |40, 41]. Spontaneous symmetry breaking can then 
be described simply by nonzero expectation values of suit- 
able bosonic fields. In a purely fermionic description the 
four fermion coupling often diverges for temperatures below 
a pseudo-critical temperature Tpc- This indicates the onset of 
local order for T < Tpc- In the partially bosonized version this 
simply translates to a vanishing "mass parameter" or gap for 
the bosons, and does not constitute an obstacle for investiga- 
tions at r < Tpc- 

Antiferromagnetic order for the two-dimensional Hubbard 
model at half filling has been successfully described in this 
framework |36, 37J. One investigates a bosonic field a which 
represents an antiferromagnetic fermion bilinear 

a(X) = i/,\X)a,Jf(X)e"^ . (3) 

We use the Matsubara formalism with Euclidean time r com- 
pactified on a torus with circumference = T ' . The Matsub- 
ara frequencies for the fermions are w = (2« + l)7rT, « E Z. 
Bosonic fields obey periodic boundary conditions such that 
cj = InnT . We use a compact notation X - (r, x), Q - (lo, q) 
and 

Z = /<"Z- IL'^t f^- 

(5(X-r) = (5(T-T')5x,x', 

m - Q') = y85„x(2^)''^^"(q - q') ■ (4) 

All components of X or Q are measured in units of the lattice 
distance a or a The discreteness of the lattice is reflected 
by the 2;r-periodicity of the momenta q. The momentum 11 in 
Eq. (O is given by 

n = (0,7r,;r). (5) 

Antiferromagnetic order is indicated by a constant nonzero 
expectation value (a(X)) = ao. The simplest description em- 
ploys a quartic effective potential for a 

, 1-9 

t/a[a] = mla+-Aaa\ (6) 

with a - a?- 12. We will see that a vanishing of corresponds 
to a diverging effective four fermion coupling, while negative 
leads to a minimum of Ua at cq + 0, and therefore in- 
dicates spontaneous symmetry breaking. In the renormaliza- 
tion group treatment, ;7i^ and ~Aa become fc-dependent running 
couplings. A situation with m^(fc) < for ksR < k < ksss, 
ml(k < fess) > describes local order in domains with linear 
scale between k^^^^ and A:^^, while no global antiferromag- 
netism is present. Macroscopic antiferromagnetic order can 
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be observed if the ^-dependent location of the potential mini- 
mum, Q'o(fc), stays nonzero as k^^ reaches the size of the exper- 
imental probe. For < a crucial ingredient for the deter- 
mination of (To is the quartic bosonic coupling In a purely 
fermionic language this corresponds to an eight fermion ver- 
tex and is quite difficult to access. 

We start with a Yukawa like ansatz for the effective average 
action 



(7) 



It describes fermion fields i// and the "antiferromagnetic boson 
field" a, with;if - (a, i//, ip*). The fermionic kinetic term 



(8) 



(9) 



involves the inverse fermion propagator 
where 



^(q) - -ju- 2f(cos qi + cos ^2) - 4f' cos qi cosq2 ■ (10) 

This is the classical inverse propagator for the Hubbard model 
with next neighbor hopping t and diagonal hopping f. The 
chemical potential is denoted by yu. We have neither included 
self-energy corrections for the fermionic propagator, whose 
significance is still under debate ||42[ 43i l44tl . nor possible 
Pomeranchuk instabilities, see e. g. i2M- In the present paper 
we restrict ourselves to a parameter and scale range where the 
spin correlations are maximal at the commensurate antiferro- 
magnetic wave vector (|5]l, e.g., we use f' = or t' /t - 0.05. 

The purely bosonic term is described by a kinetic term and 
a local effective potential Ua, cf. Eq. (|6]l. 



ra,k = ^ y a^(-Q)Pa(QMQ) + y Ua,k[a] ■ (11) 

G X 

The kinetic term Pg involves the Q-dependent part of the in- 
verse antiferromagnetic propagator and is discussed in detail 
in the appendix. The Yukawa like interaction term couples the 
bosonic field to the fermions, 

rF>.,k ^ -haJ]^\-Q)m) (12) 

e 

= -kYj ^^^+^-Q + Q>^^^ ■ ['^"'"(e)'^'A(e')] . 

K,Q,Q' 

This simplest truncation contains five ^-dependent couplings, 
namely m^, A„, h„, as well as the wave function renormaliza- 
tion Aa and the shape parameter D which parameterize PaiQ)- 
Inserting our truncation into the exact flow equation ^ 
yields a coupled system of flow equations for m^, Aa, ha and 
Pa- This is solved numerically. We start at the microscopic 
scale k - A with initial conditions 

fhl\A = U„,, ha\A = U„, ^|a = 0, P«(0|a = 0.(13) 




FIG. 1: rf-wave form factor. 



Since for this choice the effective action is quadratic in a, it is 
easy to solve for a as a functional of t//, t//*. Reinserting this 
solution into yields the well known microscopic action for 
the Hubbard model, with four fermion coupling U given by 
U = 3U,„. Equivalently, we can use 5a = Fa in the defining 
functional integral and perform the Gaussian integration over 
the "auxiliary" field a. This shows that our model is equivalent 
to the fermionic Hubbard model ll45ll . 

The equivalence between the Hubbard model and our de- 
scription of antiferromagnetic boson exchange is exact, but 
this particular form of partial bosonization is not unique. Due 
to the possibility of a Fierz reordering of the local four fermion 
interaction one could also start with a coupling both in the 
charge and the antiferromagnetic channel on initial scale or 
include even further channels 1.45,1 . Distributing the local four 
fermion action differently into the bosonic channels would al- 
ter the relation U = 3U,„- Since the different possible ver- 
sions of partial bosonization are all equivalent to the Hubbard 
model, one may question the reliability of our quantitative re- 
sults based on U,„ = U/3. In fact, mean field theory shows a 
strong dependence of the phase diagram on the choice of par- 
tial bosonization |45]. However, mean field theory neglects 
the effects of bosonic fluctuations and their inclusion substan- 
tially reduces the dependence of the results on the choice of 
bosonization [36, 46]. (Without truncations all exact partial 
bosonizations should be exactly equivalent, such that the de- 
pendence on the choice can be used as a check of the validity 
of approximations.) In this paper, we concentrate on antifer- 
romagnetic and d-wa\e fluctuations and do not include the 
charge channel. The inclusion of charge fluctuations within 
our bosonized language is in principle possible, and a first 
study 1 36] has indeed revealed that the results depend only 
weakly on the initial distribution of the four fermion interac- 
tion into the antiferromagnetic and charge density wave chan- 
nels. 

In this paper we want to understand the physics associ- 
ated to fermion bilinears whose expectation value describes 
li-wave superconductivity. We therefore have to extend the 
simplest truncation (|7]l. We introduce in addition a bosonic 
field d associated to the appropriate Cooper-pairs d 
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d(X) = -iA^(X- '-l)e^(X+'j-) - ^-ip^iX- |)e^(X+ 1). (14) 

Here, ei, 62 are the unit vectors in the plane and e = 10-2- (Eq. 
( fT4l ) is used here only as a shorthand for the definition of d 
in momentum space, Eq. ( l34b . since the Grassmann variables 
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are not located at the lattice sites. A definition of d in position 
space and the appropriate Fourier transformation can be found 
e.g. in a.) 

For the purely bosonic part we add in our truncation solely 
a mass term of the cZ-boson 

r^^'w^X^'^fiMG). (15) 

G 

The £/-field couples to the fermions by a Yukawa term 
rFd,k = -hdY,(d*(X)d(X) + d(X)d*(X)) (16) 

X 

= J](d*(Q)d(Q) + d{Q)d\Q)) 
Q 

K,Q,Q' 

x[d*{K)[V{Q)em')] - d{K)[ilj\Q)er{Q')]), 
with Yukawa coupling lid- The cZ-wave form factor 

/rf(q) S COS y - COS y , (17) 

is kept fixed and shown in Fig. [T] We note the characteristic 
change of the sign under rotations of 90°. In our approxima- 
tion, only the momentum independent coupling hd depends on 
the scale k. The extended truncation has two further running 
couplings, and /zj. 

The initial values in the d-wave channel are 

mlk^l, hdlA^O. (18) 

At the microscopic scale the li-boson therefore decouples 
from the fermions and the a-boson, such that the microscopic 
action is not modified. Our model remains exactly equiva- 
lent to the Hubbard model. For k < A, however, a non-zero 
Yukawa coupling hd is generated, as shown in Fig. |2] "In- 
tegrating out" the c/-boson by solving the field equations as a 
functional of i//, and reinserting into Ft, yields now an eff'ec- 
tive four fermion interaction in the li-wave pairing channel 

Even though absent microscopically, this interaction is gen- 
erated during the flow by the coupling to the antiferromag- 
netic channel, as derived before in a purely fermionic lan- 
guage illlHlllHIlSIIlllillllMElt]- The eff^ective 
coupling /l^ is shown graphically in Fig. |2l 

III. FLOW EQUATIONS 

In addition to the truncation of the effective average ac- 
tion, the regulator functions have to be defined. Our 
choice for the fermionic regulator is inspired by the fact 
that at nonzero temperature the inverse fermionic propagator 



0.5 



0.4- 



0.1 - / 

0.0 r ' ' ^-r'^-^ - 

-2 2 4 6 8 10 

-In^/f 

FIG. 2: Generation of d-wave coupling. The solid line shows the 
flow of the d-wave channel fermionic coupling /l^ for T ft = 0.13, 
fi/t = -0.1 and t' = 0. The dashed line shows the flow of the Yukawa 
coupling in the rf-wave channel /ij, as given by rebosonization for the 
same T, ji, f and initial conditions dlSt . 

Pf(Q) - 2ni{nf + j)T + has no zero eigenvalue. The tem- 
perature itself acts as a regulator. We put this into use by the 
regulator function lf36ll 

< (0 = ic^iy - 1) = 27r/(nf + i)(r, - T) , (20) 

with 

= T'^ + k'^, (21) 

-r /; /-r \3 f 1 if » T 

d,n - (k/n) ^ I ^^^y,^3 ^^^^ . 

For k > T the cutoff Rj^ in the inverse fermion propagator 
suppresses the contribution of all fluctuations with momenta 
Iq - p < (nk)^, even for T = 0. It becomes ineffective for 
k where no cutoff is needed anyhow. Basically, the tem- 
perature T is replaced by the scale dependent "temperature" 
Tk - we cool the fermions down to the temperature of interest 
during the flow. For k ^ oo the cutoff diverges for all T such 
that fermion fluctuations are completely suppressed and we 
start with the same initial conditions for all T . For A: — > the 
cutoff function vanishes for all T and one recovers the quan- 
tum effective action. For finite A:"', corresponding to a finite 
macroscopic size of the probe or a finite experimental wave 
length, the temperature dependence of the cutoff is small for 
A: <K r. In principal, different functional forms of Tk{k) can 
be used for a test of the robustness of our truncation. 

For bosons the situation is different. Here, long range 
bosonic modes may cause infrared problems which cannot be 
regularized by the same type of regulator function we use for 
the fermions. In particular, near a second order phase transi- 
tion the bosonic correlation length diverges, which is the same 
as a vanishing boson mass term. Our regulator is devised 
in order to cut off the long range bosonic fluctuations. We 
specify the regularization of the a-bosons in detail in the ap- 
pendix. No regularization of the li-bosons is needed for the 
investigations in the present paper 
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FIG. 3: One-loop corrections (box diagrams) to the four fermion 
interaction. Solid lines represent fermions, wiggly lines a-bosons. 



The flow equations for the couplings follow from projection 
of the flow equation (|2| onto the corresponding monomial of 
fields. For the mass parameter in the effective potential one 
derives Il36ll 



1 



Q pMQ)pMQ + n) 

5- -y- 1 



(22) 



where 



P'a(Q) = Pa{Q) + KiQ) , pUQ) = Pf(Q) + RliQ) ■ (23) 

The derivative dk acts only onto the explicit A;-dependence in- 
troduced by the regulator functions, not on the couplings. For 
the quartic coupling the flow equation reads 



«Zi *(p^^(0p^^(e + n))2 



(24) 
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1 



The flow of the kinetic term PaiQ) is discussed in the ap- 
pendix. 

The flow equation for the Yukawa coupling consists of a di- 
rect contribution /3j^ and a contribution /3j^ which comes from 
rebosonization of regenerated fermionic couplings H 



The first contribution reads 



1 



pWPf^Q + n)(nxe) + ml) 



(25) 



(26) 



The second contribution is discussed in the foHowing section, 
together with the flow equation for hj. For the mass parameter 



for the li-boson we find 



/i(2q) 



(27) 



IV. GENERATION OF d-WAVE COUPLING 

In this section we discuss how a coupling in the li-wave 
channel is generated by the exchange of antiferromagnetic 



bosons. This mechanism was suggested e.g. by lUlllll]. At the 
microscopic scale A we have a local repulsive four fermion 
interaction, described in our picture by the exchange of a- 
bosons. No other channel is present. However, even in the 
case of purely repulsive forces between the electrons the sys - 
tem may become unstable against Cooper-pair formation IHTIl . 
On lower momentum scales we indeed observe a generation 
of an attractive coupling in a Cooper-pair channel which turns 
out to have li-wave symmetry. This may become critical in an 
appropriate parameter region. 

We show in Fig. [3]the box diagrams exchanging a-bosons. 
They generate an effective momentum dependent four fermion 



vertex F 



<4) 
F 



Q\,-,Qi 



Xra(Ql)MQ2Wy(Q3)MQ4) ■ (28) 



This irreducible vertex has to be added to the effective interac- 
tion arising from the exchange of bosons if we account for the 
total effective action between the electrons. It receives further 
contributions beyond the box diagrams since it has to include 
all contributions to the total four fermion vertex which are not 
accounted for by the antiferromagnetic boson exchange in a 
given truncation. In particular, the fluctuations of the antifer- 
romagnetic bosons induce a momentum dependent piece in 
the effective Yukawa coupling between the electrons and the 
a-bosons. This momentum dependent part is not reproduced 
by our truncation where a constant h„ is evaluated for a partic- 
ular choice of external momenta. We have therefore added this 
piece to Our computation of the sum of r4 and the anti- 
ferromagnetic boson exchange agrees with a purely fermionic 
computation in one loop order, such that the results are iden- 
tical as long as U/t remains small. 

The li-wave channel interaction is part ofr^^^Qf). Since the 
coupling r^'^\Qi) contains contributions from various chan- 
nels we will next project onto its contributions from the d- 
wave and the antiferromagnetic channels. For spin rotation 
invariant systems the spin structure of the four fermion inter- 
action has the general form 



^'^^(Gi, 22, 23, 24) = r^":(2i,22,23,24)5„y;/,^ (29) 



<4), 



+4*)(2i,22,23,24)Wtf, 

where S ayfie - Sa/iSyg - 5yp6a6 and Tay^ps - 5ap5y6 + Sypdas 
project onto spin singlet and spin triplet states in the fermion 
pair channels, respectively. Here, we are interested in spin 
singlet states and we choose Q\ - -Qi = L - (nT, 1) and 
Qi - -Qa = L! - {ttT, V). For the projection onto the d-wa\e 
Cooper-pair channel we use the fact that the d-wa\e coupling 
changes its sign under rotation of 90°. We therefore define a 
momentum dependent d-wave channel coupling Ap{l, I') by 



4(1,1') = ~{rfl(L,L',-L,-L') 

-I<^l(R(L),L',-R(L),-L')], 



(30) 



where R{L) denotes a rotation of the spatial components 1 of 
L by 90°. For the definition ( l30b we have kept L' fixed and 



I -n ^ Z -n 

(a) T/t = 0.5 , ^i/t = -0.1 (b) T/t = 0.5 , ^J^.|t = -0.1 




FIG. 4: Momentum dependence of rf-wave coupling. In (a), (c) 
we show the normaUzed momentum dependence of the rhs of the 
flow equation for the fermionic coupling, ^((9j/l^(l, r)/5t/i^). In 
(b), (d) we display the residual coupling after the subtraction of the 
bosonized part, - ,/rf(21)/rf(21')) (cf. Eq. (O). We 

use /' = in all plots. Notice that these plots are independent of U . 

subtracted the same contribution after a rotation of 90° of the 
space-like components of L. For instance, the antiferromag- 
netic channel and the i-wave channel are subtracted in this 
way. 

We compute dkT^P and 1') by taking the fourth func- 

tional derivative of Eq. (|2]l with respect to the fermion fields. 
Within our truncation this corresponds to the box diagrams in 
Fig. [3] We find that 1') is well approximated by the 

flow of a simple cZ-wave channel coupling /l^ of the form 

4(l,l') = /rf(21)/rf(21')4. (31) 

In Figs. |4^,c we show the flow of the Ihs of Eq. ( ISTT i as a 
function of 1 for 1' = (0, tt) (for k - K and with appropriate 
overall normalization). There is a good qualitative agreement 
with the d-wa\e form factor displayed in Fig. [1] We determine 
the flow of 4 by choosing 1 = (tt, 0), 1' = (0, n), i.e. 5^/1^ = 
5,4((;r,0),(0,;r)). 

Next, we show how this fermionic quartic coupling can be 
bosonized on all scales k. This means that we trade the four 
fermion coupling in favour of a Yukawa-coupling of the 
li-bosons. After this transformation, the li-boson exchange in- 
duces the effective interaction in the li-wave pairing channel. 
We achieve this reformulation by a variable transformation in 
the functional flow equation The flow equation (O de- 
scribes the scale dependence of at fixed fields x- We will 
now switch to new fc-dependent variables Xk and derive new 
flow equations which describe the ^-dependence of at fixed 
Xk- For a suitable choice of Xk this will realize the desired 
bosonization |34, 35]. 

More precisely, we choose scale dependent bosonic fields 
B/t - Bk[^,iff* ,B; k] = (ak,dk,d*^) and perform the variable 



FIG. 5: Effective four fermion interaction in the rf-wave pairing chan- 
nel mediated by d-boson exchange. 



transformation in the exact flow equation 

^r,[^,^%B,]|3^ = 5ir,[«A,iA*,B]|3^3^ (32) 
+ Z(^r,[^,^-,B,])a,a,(0 

+ Z(^r,[^,r,B,])4«'*(0 

+ Z(^r,[^,r,B,])5K(e). 

This equation is our starting point for the construction of the 
"perfect bosons" [34J on all scales. The first term on the rhs 
is given by Eq. (O. 

The scale dependent boson fields are defined by adding suit- 
able fermion bilinears 

am ^ adQ) - atm) , (33) 
dk(Q) ^ ddQ) - 4d(Q) ■ 

Here the fermion bilinears 

a(e) = 2f'"(P)(TiA(P-n + 0, (34) 

p 

d(Q) = i2/d(2p-q).AV)e<A(G-^), 

^ p 

correspond to the Fourier transforms of Eqs. (O, (fl4l i. In 
Eq. (l33T l. we denote by and the original microscopic 
fields appearing in the microscopic action Fa. These are the 
A:-independent fields which are kept fixed in Eq. (|2]l or in the 
first term on the rhs of Eq. (l32b . The fermion fields are kept 
scale independent and we infer 

dkak(Q) = -dkaim), (35) 
dkdf>(Q) = -duaid^'KO). 



The variable transformation in Eq. ( [32] l induces new 
Yukawa and four fermion couplings, given in our truncation 
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by 



e 

Q 

+ Y,idkal)hMQ)-K-Q) 

Q 

+2Yj.dkai)hdd\Q)d{Q) . (36) 



Here the couplings hg, hd etc. correspond to suitable func- 
tional derivatives of at fixed Bj. 

For a suitable choice of a^, we can use the four fermion in- 
teraction in the d-wave pairing channel generated by the vari- 
able transformation, in order to cancel a corresponding piece 
generated by the diagrams of Fig. [3] in ^tr;; . More 



precisely, we fix dko'l by the requirement 



O^dkAi+lhddua'l, 



(37) 



where the first term is computed from Eqs. (|2]i, (l28Tl-(l3ni. The 
coupling /l^lg (1,1') is defined similar to Eq. ( |30] |. However, 
it is now given by the functional derivative of F^: at fixed B^:. 
This means that its flow is obtained from the flow of /l^(l, 1') 
(at fixed B) by subtracting the "bosonized piece" 

5,4(1, 1')|b^_ = 5*4(1' 1') - dk4fd{2\)fd{2\') . (38) 

In particular, if we define in analogy to Eq. (ISTT i. one 

has <5i4lBi " <9i4(l, 1')|b, i" Figs- Ht-d (in a 

suitable normalization where 5^4 divided out). Compari- 
son with Figs. |4^,c shows that the major part of 5^4(1' ' ) 
indeed be absorbed by the bosonization. The "residual part" 
3a:4(1' 1')|b contains higher harmonics - the dominant piece 
in Figs. |4j5,d has a periodicity in momenta with period n, in- 
stead of the leading piece with periodicity 2k in Figs. |4^,c. 
In principle, this residual interaction can be partially absorbed 
by a more complex fc-dependent field in Eq. ( |33] l. We will not 
do this here and simply neglect in our truncation the residual 
interaction in (9<:4(1' 1')|b ■ '-^'^ ^^^^ ^ comparison of 
Figs. |4j) and|4}l, the accuracy of this approximation gets bet- 
ter for larger |yu/r|. After the variable transformation, we can 
work with a vanishing four fermion coupling in the li-wave 
pairing channel. 

From Eq. ( |32] | we can now infer the additional contribu- 
tions to the flow of the Yukawa coupling lid- It accounts for 
the fact that after the variable change a leading contribution 
to the four fermion interaction is mediated by the exchange of 
bosons and no longer contained in /If . One obtains a modified 
flow equation for the Yukawa coupling for the d-boson 



dkk 



knd 



dkh 



'knd , 



'^OkOk 



Thereby inserting ^kC^^^ from Eq. (l37l l yields 

mldkXi 



dkhd = dkhdL - ^, 
2h 



(39) 



(40) 




FIG. 6: Scale dependence of ml, in^ and A„ for temperature above 
Tp, {Tit = 0.12, solid lines) and below Tp, {Tjt = 0.09, dashed 
lines). We use ult = -0.1, f = 0, Ujt = 3. 



for the li-boson. In our truncation one has dkhd\^ - 0. The 
numerical solution to this flow equation is shown in Fig. |2] 
We can reconstruct the effective four fermion interaction in the 
t/-wave pairing channel by computing the tree diagram from 
the (i-boson exchange (cf. Fig. |5]l 

4eff(ei' 22, a, Qa) = -%fd{<\i - q3)/d(q2 - q4) (41) 

xS(Qi -Q2 + Q3- Q4) ■ 

Similar as for the c/-boson, the Yukawa coupling of the 
a-boson receives an additional contribution due to the re- 
bosonization of the four fermion interaction in the antiferro- 
magnetic channel ll36ll . 



dkha = dkhJ + mldkal . 



(42) 



We define the projection onto the four fermion interaction in 
the antiferromagnetic channel by 



4 = M%i(o,n,o,-n), 



and determine dka'l by 



= 5,4 + hAal . 



This yields 



— — 

dkha = dkha\^ - -fdkA"^ +Pl 

tin ' 



.rh 



(43) 



(44) 



(45) 



for the flow of the Yukawa coupling of the a-boson. The 
first contribution was given in Eq. ( l26l l. The second contri- 
bution comes from rebosonization and we find in agreement 
with S 



1 



X 



{Pi{Q) + ml){Pi(n-Q) + ml) 

1 



(46) 
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FIG. 7: ln(-4//) for yt = is plotted as a function of T/f for T 
Tpc = 0.1004f, yilt = -0.1, /' = 0. 



FIG. 8: ln(-^° //) for yt = is plotted as a function of T/t for T 
Tp, = 0.1004r, yu// = -0.1, f = 0. 



We have solved the flow equations for the couplings ha, hd, 
ml, iiij, Aa and D numerically. In Fig. |6] we show the k- 
dependence of m^, ilij and Ag. For T > Tpc both and 
remain positive for all k. We show the temperature depen- 
dence of the effective four fermion couplings in the li-wave 
and the antiferromagnetic channel, /l^ and A"p, in Figs. |2]and 
[8] Both couplings are evaluated for ^ = 0. 



V. PSEUDOCRITICAL TEMPERATURE 

For T — > Tpc the "mass term" reaches zero at a nonzero 
value of k, and therefore A"^ diverges. This feature extends to 
all T < Tpc and signals the onset of local antiferromagnetic 
order in domains of linear size A; For smaller k, the flow 
has to be continued in the "spontaneously broken regime" |l 36l 
[37I1 . where the minimum of the effective potential occurs for 
nonzero a. In the present paper we do not continue the flow 
beyond the scale where ml vanishes, but rather stop once ml 
reaches zero. The corresponding pseudocritical temperature 
Tpc is shown in Fig. |9]as a function of the chemical potential. 
We also show the dependence of Tpc on the Hubbard coupling 
U in Fig. [To] For U/t = 3 and ju = f' = we find for Tpc a 
value that is about 30% smaller than in |36]. This is due to the 
different ansatz for the propagator of the a-boson and may be 
taken as an indication for the size of the error. In the present 
work we restrict ourselves to Hubbard couplings U/t < 3 and 
to temperature T/t > 0.02. In order to extend the lines in Fig. 
[To] to smaller T one has to improve our simple truncation of 
the frequency dependence of the propagator of the a-boson. 
This issue will be addressed in a forthcoming work. 

Another interesting quantity is the relative strength of the 
effective four fermion interaction in the li-wave pairing and 
the antiferromagnetic channel. We plot the ratio 



iP.ihl 

a a 

hlml 



(47) 



in Fig. [TT]for different values of //. In these plots the ratio 
R first increases due to the generation of a non-vanishing lid, 
and subsequently decreases since ihl approaches zero or small 



values. We see the tendency that the maximum of R increases 
with increasing //. 

Our present truncation does not remain valid for low k if 
\^/t\ becomes large. The reason is that Ag or Aa reach zero 
before either ml or m^^ vanishes. This means that our approx- 
imation for the momentum dependence in PaiQ), or for the 
polynomial approximation for Ua, becomes insufficient. An 
extended truncation can cure these problems, which will be 
addressed in a future publication. In our figures the break- 
down of the approximations can be seen by the end of the 
curves in Fig. [9] and Fig. [TT] The insufficiency of the trunca- 
tion is the reason why we cannot access the values of |ju| and 
k where R becomes large in the present paper. 



VI. CONCLUSIONS 

We have shown how the renormalization flow turns the mi- 
croscopic Hubbard model into an effective boson exchange 
model at intermediate length scales k^^. This "mesoscopic" 
model is characterized by the exchange of collective electron 
or electron-hole pairs, which are described by bosonic fields. 
The most important channels are the exchange of electron- 
hole pairs of the antiferromagnetic type - the a-boson - and 
li-wave electron pairs - the J-boson. The coupling of these 
bosons to the electrons is described by Yukawa type cou- 
plings ha, hd, which multiply the form factors appropriate 
for the respective channels. The effective model is char- 
acterized by the effective action, cf. Eqs. (O, dTSl l, ( [T6b . 
r = Tf -H r„ -H Tfa +rd + rfd. Besides the Yukawa couplings, 
important ingredients are the propagator and the effective po- 
tential for the a-bosons and the li-bosons. 

The effective boson exchange model has several coupling 
constants, as ha, hd, ihl, well as bosonic self inter- 

actions and couplings parameterizing the bosonic propaga- 
tors (cf. appendix). These couplings depend on the mo- 
mentum scale k at which the effective model is considered. 
We may consider the effective boson exchange model as a 
coarse grained version of the microscopic Hubbard model, 
with coarse graining length A; Its couplings can therefore 
be related to the couplings f, t' and U of the Hubbard model. 
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FIG. 9: Pseudo-critical temperature Tpdt as a function of yu/f for 
U/t = 3 and /' = (solid line) and t' = -0.05/ (dashed line). 



and we have done so by the use of the functional renormahza- 
tion flow. This flow shows how an effective coupUng in the 
c/-wave pairing channel is generated by the fluctuations of the 
a-bosons. 

One may also consider the effective boson exchange model 
with arbitrary parameters. This can then be interpreted as 
a generalization of the Hubbard model. The particular case 
ha = 0, where the a-bosons play no role, has already been 
discussed with the help of functional flow equations in Ref. 
ifssll . For low temperatures one finds superfluidity. The in- 
vestigation of pSJ] has covered the critical behavior at the 
phase transition and the low temperature "ordered phase". 
The phase transition was found to be in the universality class 
of the Kosterlitz-Thouless phase transition. Reinterpreted in 
the context of our investigation of the Hubbard model, the 
model 1 38] becomes a valid approximation if at some scale 
k the antiferromagnetic channel coupling /l^ = -hi/ 2ml be- 
comes much smaller than the li-wave coupling /l^ - -h^/fhj. 
The corresponding large values of R are reached, however, 
only in a region for /i/f where our present approximation for 
the propagator of the antiferromagnetic bosons becomes in- 
sufficient. 

As we have mentioned in Sect. V, our present investigation 
is limited to the range of k where the propagator for the 
a-boson is well described by the approximation discussed 
in the appendix. In particular, this requires Aa > 0. We 
implicitly assume that commensurate antiferromagnetic 
fluctuations are dominant as compared to incommensurate 
antiferromagnetism or ferromagnetism. The latter play 
an irnportant role in other parameter ranges of this model 
S S m H. Wifliin file present truncation the 
positivity of Ag is not realized for all jj. ,t ,t' ,U ,T if k 
becomes small. Also the simple dependence of Pa on the 
Matsubara frequency does not remain appropriate for T — > 0. 
We will discuss a more general form of Pa, as well as a 
momentum dependent propagator for the li-boson, in a 
forthcoming work. We are confident that a suitable trunca- 
tion of the functional renormalization group equations will 
give access to the whole phase diagram of the Hubbard model. 

Acknowledgments: We thank S. Diehl, J. M. Pawlowski, 
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FIG. 10: Pseudo-critical temperature Tpc/t as a function of U/t for 
^ = 0,f' = (solid line) and n/t = -0.12,/'// = -0.05 (dashed 
line). 
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APPENDIX: BOSONIC PROPAGATOR 

In order to gain some first information on the general shape 
of the momentum dependent piece in the inverse propagator 
for the antiferromagnetic boson, we compute the contribution 
from the fermionic loop 

^^"(2) = Z p,,Q,pln)PAP) " ■ 

The general features of this one loop contribution will be used 
in order to motivate a suitable truncation. Performing the Mat- 
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FIG. 1 1 : Relative strength of li-wave and antiferromagnetic coupling. 
The ratio R, defined in Eq. l i47! > is plotted for four parameter choices: 
1.) T/t = 0.1235,yL( = 0, 2.) T/t = 0.115,//// = -0.09, 3.) T/t = 
0.08,//// = -0.15 and 4.) T/t = 0.05,//// = -0.3. We use U/t = 3, 
/' =0. 
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FIG. 12: Mean field calculation of the frequency dependence of the 
bosonic propagator at zero spatial momenta, m^U + AG^^(2mnT, q = 
0) for U/t = 3, T/t = 0.25, f' = and = 0. 



subara sum one finds 



2 (2n 



d^p tanh(§) - tanh(i^ + f ) 



q+p+;r 



— ia> 



(A 2) 



In Fig. [T2] we plot the frequency dependence at vanishing 
spatial momenta while the dependence on spatial momenta at 
vanishing frequency is plotted in Fig. [T3](a). 

We approximate the inverse bosonic propagator at zero fre- 
quency by 

rQ]2r)2 

P„,,(0, q) = AG,-)(0 + m^U - m% = A,-^^ ^ (A 3) 

where [q]^ is defined as [q]^ - l\ + I2 ^'^^ 1' ^ [~^'^] ™d 
continued periodically otherwise. It is convenient to define 
the gradient coefficient Ag and the shape coefficient D by 



and 



A. = i^P«(0,/,Ot„ 



D^^-^iPa(0,7:,7T)-P,(0,0,0))- 



(A 4) 



(A 5) 



Comparison between Figs. [T3h and [T3b shows that this 
choice reproduces the mean field result rather well. (Note that 
Pfl(0, n, 71) computed from Eqs. dA 3l l - jA 5l l deviates slightly 
from the mean field result for this momentum choice.) The 
lowest order perturbative contribution from the fermions reads 



52 



df P'p(Q)P''p(K +Q + n) 



(A 6) 



For frequencies different from zero we set for all q 



Pa{cL> * 0,q) ^ ml\A - milk . 



(A 7) 



In practice, this means that we include the fluctuation effects 
only for the zero Matsubara frequency mode. In the temper- 
ature range of interest in this work, this is a well justified ap- 
proximation, cf. Fig. [12] On the other hand, the truncation 
(lA 71 ) is not suited for very low T and therefore limits the tem- 
perature range for our investigation in the present paper 

Within the functional renormalization group approach, we 
describe the scale dependence of the inverse bosonic propaga- 
tor by flow equations for the parameters Aa and D. The flow 
equation for the gradient coefficient is 



^ " 'dfpX^{Q)pX^{K 



e + n) 



(A 8) 



/=o 



where K - (0, 1, 0). At the initial scale the gradient coeffi- 
cient vanishes Aa\h - 0. The boson is not dynamic on this 
scale and can be regarded as an auxiliary field. The difference 
Pa{Q,n,ji) - fo(0, 0, 0) in Eq. dA 5l l is also scale dependent. 
From this, and Eqs. ( lA 51 ). ( lA 81 ). one infers the flow equation 
for D. Again, D vanishes at the initial scale A. The result 
of our approximation dA 31 ), as calculated from the flow equa- 
tion, is displayed in Fig. [T3]b. Comparison with the one loop 
result in Fig. [T3]a shows satisfactory agreement. 

The infrared cutoff 7?^ for the a- boson i s _adapted to Eq. 
dA 31 ). We use an "optimized cutoff" 



D2 + [q] 



KiQ) = Aa ■ {e - 



D2 + [q] 



r). (A 9) 
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FIG. 13: The mean field approximation for the spatial momen- 
tum dependence of the bosonic propagator at zero frequency m^U + 
AG;'(0,q) is shown in (a). We use U/t = 3, yu = 0, f' = and 
T/t = 0.25. In (b) the parameterization in terms of Eq. dA 3t is 
displayed for a renormalization scale ^ = 0. 
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